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ABSTRACT 

The atmospheric properties above three sites (Dome C, Dome A and the South Pole) 
on the Internal Antarctic Plateau are investigated for astronomical applications us- 
ing the monthly median of the analyses from ECMWF (the European Centre for 
Medium-Range Weather Forecasts). Radiosoundings extended on a yearly time scale 
at the South Pole and Dome C are used to quantify the reliability of the ECMWF 
analyses in the free atmosphere as well as in the boundary and surface layers, and to 
characterize the median wind speed in the first 100 m above the two sites. Thermody- 
namic instability properties in the free atmosphere above the three sites are quantified 
with monthly median values of the Richardson number. We find that the probability to 
trigger thermodynamic instabilities above 100 m is smaller on the Internal Antarctic 
Plateau than on mid- latitude sites. In spite of the generally more stable atmospheric 
conditions of the Antarctic sites compared to mid-latitude sites, Dome C shows worse 
thermodynamic instability conditions than those predicted above the South Pole and 
Dome A above 100 m. A rank of the Antarctic sites done with respect to the strength 
of the wind speed in the free atmosphere (ECMWF analyses) as well as the wind shear 
in the surface layer (radiosoundings) is presented. 
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1 INTRODUCTION 

The summits of the Antarctic Plateau are of particular in- 
terest for ground-based astronomy because the optical tur- 
bulence appears to be confined to a narrow layer close 
to the ground, (|Marks et al.l Il999l ; lLawrence et al.1 12004 



lAgabi et al.ll2006h . Above this layer the atm osphere is excep- 
tional ly clea r and the turbulenc e weak (0.27' lLawrence et al.l 
i|2004h . 0.36' lAgabi et al l i|2006l ) at Dome C). Measurements 
have shown that the height of the turbulent layer above 
the sum mits is much lower than above other sites on the 
Plateau. (Lawrence et al.1 12004 iMarks et all 1 19961 : iMarksl 
120021 ; IMarks et al.1 1 19991 ). More precisely, the height of this 
layer is much lar ger above the South Po le (220 m (|Markj 
I2OO2I ) or 270 m jTravouillon et all I2003T0. which l i es on a 
slope, than above Dome C (36 ±10 m lAgabi et~afl (|2006t )1. 
Above Dome A the turbulence is expected to be even weaker. 

The surface winds of Antarctica are among the prin- 
cipal sources of turbulence near the ground. The dominant 
source of the surface winds is t he sloping terrain an d the 
radiative cooling of the surface (|Schwerdtfegerl 11984) . The 
radiative cooling produces a temperature inversion that to- 
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gether with the sloping terrain cause a horizontal temper- 
ature gradient. This triggers a surface wind alongside the 
slope of the terrain. Strong wind shears can therefore occur 
in the boundary between the surface winds and the winds in 
the free atmosphere, which in general are geostrophic. This 
is the main source of instability under conditions of extreme 
stability, as is the case of the Antarctic atmosphere in win- 
ter. Above the summits of the Internal Antarctic Plateau 
the surface winds should be much weaker than elsewhere on 
the Plateau due to a lack of the principal cause triggering 
them: a sloping terrain. 

These elements justify t he enthusias t ic int e rest of 
astro nomers for this site l)Storev et alJ 120031 : iFossatl 
120051 ) Studies of the atmospheric properties of these regions 
are fundamental for applications to ground-based astron- 
omy. However we need a better quantification of these char- 
acteristics, using instruments and measurements as well as 
models and simulations, in order t o fill the gaps of unce r- 
tainties or doubts that still remain (|Masciadri et al.ll2006l ) , 
extending our attention to a comparative analysis of differ- 
ent sites of the Internal Antarctic Plateau. 

It is important to produce statistical estimate of me- 
teorological parameters at all heights from the ground to 
verify if atmospheric conditions are always advantageous 
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for astr onomical applications. Ind eed, it has recently been 
shown (Geissler & Masciadri 2006), using European Centre 
for Medium-Range Weather Forecasts (ECMWF) analyses, 
that in winter the wind speed grows monotonically above 
~10 km a.s.l. (above sea level), achieving median values of 
the order of ~30 m s _1 at 25 km a.s.l. At this height a vari- 
ation of ~20 m s~ has been estimated between summer 
and winter. Such a strong wind might trigger an important 
decrease of the wavefront coherence t ime, even in presence 
of a w eak turbulence (see discussion in Geis sler fc Masciadril 
(2006)) and, as a consequence, the potentiality of these sites 
might vanish. It should be therefore interesting to better 
quantify the median wind speed profile on other sites (with 
astronomical interest) of the Internal Antarctic Plateau or 
to retrieve some general indication of the wind speed in the 
free atmosphere above the Internal Antarctic Plateau. 

Besides that, the employment of ECMWF analyses for 
char acterization of the sur f ace la yer requires a deeper anal- 
ysis. IGeissler fc Masciadril (|2006l ) concentrated their study 
at heights greater than 30 m thus excluding the surface 
contribution, assuming that the ECMWF analyses are not 
optimized for the atmospheric flow near the surface. More 
recently, studies l|Sadibekova et al] 120061 ) appeared claim- 
ing that the ECMWF analyses can be used to quantify and 
characterize the atmospheric properties with a good level of 
accuracy down to the surface level. In spite of our conviction 
that this conclusion is the resu lt of a partial analysi s (only 
summer data) we admit that in lGeissler fc Masciadril (120061 ) 
the authors assumed (and they did not proved) the limits of 
the ECMWF for the surface layer. We therefore think that is 
time to provide a dedicated analysis on this subject to know 
the limit within which we can achieve reliable estimates with 
General Circulation Models (GCM) i.e. with ECMWF anal- 
yses and to identify the domains in which one is forced to 
employ mesoscale models. The latter are in principle more 
suitable to better resolve phenomena happening at smaller 
spatial and temporal scales. 

The usefulness of mesoscale models depends on the lim- 
itations imposed by the products of the General Circulation 
Models (that means the ECMWF analyses). It is obvious 
that, the usefulness of mesoscale models would not be justi- 
fied for wind speed if ECWMF products could provide an- 
swers with a sufficient good accuracy. 

Our group is involved in a long term study made with 
meso-scale models for the simulation of the optical turbu- 
lence on the whole tropos phere and low stratosp here using 
the technique described in [Masciadri et al.l (|2004r ) above the 
Internal Antarctic Plateau_J. This study is therefore prope- 
deutic to researches done with such a typology of models. 

It is therefore fundamental to provide a clear picture of 
the limitations of the ECMWF products and at the same 
time, to try to retrieve the maximum number of useful in- 
formation we can get from such a kind of products. 

In this paper we try a first attempt to quantify, above 
the three sites with some astronomic interests (the South 
Pole, Dome C and Dome A), the differences of some criti- 
cal meteorologic parameters that are directly, or indirectly, 
related to the characteristics of atmospheric turbulence. We 

Som e attempts have been done in the past llSwain fc Galled 
2006) even if with different scientific goals. 



Table 1. The geographic coordinates of the sites and the closest 
grid points from which the ECMWF analyses are extracted. 



Site 


Lat. 


Long. 


Dome A* 


80°22'00"S 


77°21'11"E 




80°30'00"S 


77° 30' 00" E 


Dome C ** 


75°06'04"S 


123°20'48"E 




75°00'00"S 


123°30'00"E 


South Pole 


90°00'00"S 


0°00'00"E 




90°00'00"S 


0°00'00"E 



* Measured with GPS by Dr X. Cui (private communication). 
"Measured with GPS by Prof. J. Storey (private communica- 
tion). 

expressly select two sites (South Pole and Dome C) for which 
measurements are available and one site (Dome A) for which 
no measurements are available. The reasons for this choice 
is explained later on (paper's scientific goals). 

We use data from the MARS catalogue of the ECMWF 
and radiosoundings from the South Pol^l and Dome 

cB 

Analyses data are extracted from the three grid points 
that are nearest to the sites of interest, i.e. Dome A, 
Dome C and the South Pole. The coordinates of the sites 
are given in Table [1] We extracted analyses data from 
MARS at 00:00 UTC for the whole year of 2005 to assure 
a complete statistical sample covering all seasons. A more 
detailed description of the analyses data set is given by 
IGeissler fc Masciadril (|2006m . 

The scientific goals of this paper are: 

(1) To carry out a detailed comparison of radiosound- 
ings/analyses of the wind speed and the potential temper- 
ature (the main critical parameters defining the stability of 
the atmosphere) near the surface (the first 150 m) for winter 
as well as summer above Dome C and the South Pole. This 
will permit us to quantify the uncertainty between measure- 
ments and ECMWF analyses in this vertical slab. The idea is 
therefore to define the conditions in which the ECMWF can 
be used to characterize, with a good level of reliability and 
accuracy, some atmospheric parameters and to use this tool 
to characterize a site for which no measurements are avail- 
able. This is of course the interest for a model. Depending 
on the results of this analysis we will perform comparisons 
of meteorologic parameters above the South Pole, Dome C 
and Dome A using ECMWF analyses (Section 2). 

(2) Using radiosoundings we will estimate the statistic 
median values of the wind speed in the first tens of meters 
above the South Pole and Dome C extended from April to 
November. We will therefore be able to quantify which site 
shows the better characteristics for nightly astronomical ap- 
plications (Section 3). 

(3) We extend the study developed by 
IGeissler fc Masciadril ([2006) above Dome C for the 
ECMWF wind speed also to the South Pole and Dome 
A, both located on the Internal Antarctic Plateau but at 
different latitude and longitude. In this way, we intend to 
quantify which site is the best for astronomical applications. 

2 ftp:/ /amrc. ssec.wisc.edu/pub/southpole/radiosonde 

3 http://www.climantartide.it 
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Results of this analysis are fundamental to confirm or see 
in the right perspective the potentialities of Dome C. 

(4) We extend the analys i s of t he Richardson number 
done bv lGeissler fc Masciadril l|200rj ) for Dome C for heights 
above 30 m to the three sites (South Pole, Dome C and 
Dome A) in order to quantify the regions and the periods 
that are less favorable for the triggering of optical turbu- 
lence and to identify the site with the best characteristics 
(Section 5). This result should represent the first estimate of 
potentialities of Dome A for astronomical applications and 
this should mean that we are able to provide some reliable 
results and conclusions even before some measurements are 
done on that site. This study has therefore a double inter- 
est. Firstly the intrinsic result itself. Secondly this analysis 
should open the path to a different approach for a fast and 
reliable classifications of potential astronomical sites. 



2 ECMWF ANALYSES VERSUS 
RADIOSOUNDINGS 

In this section the ECMWF analyses are compared to ra- 
diosoundings from Dome C and the South Pole in order 
to investigate the reliability of the ECMWF data over an 
antarctic site. The radiosoundings and the ECMWF anal- 
yses used for this comparison refer to the year 2006. This 
year was chosen because of the richer sample of available 
radiosoundings. The reason why the comparison of analy- 
ses from different sites, discussed later on in this paper, are 
from 2005 is because we wish to investigate a homogeneous 
data set. In February 2006 the number of vertical levels in 
the ECMWF model changed from 60 to 91. 

When studying the difference between radiosound- 
ings and analyses particular interest is paid to the first 
tens of meters since this is where we expect the largest 
turbulent activity and th is range was not studied by 
iGeissler fc Masciadril l|2006T ). The median of the difference 
in wind speed and temperature between ECMWF analyses 
and radiosoundings, for summer (December, January and 
February) and winter (June, July and August) are shown in 
Fig. [T] A total number of 73 nights has been used in summer 
and 75 in winter. In the free atmosphere we observe that the 
radiosoundings do not reach as high altitudes in winter as 
they do in summer. A similar effect ha s been observed above 
the S outh Pole in our previous study (IGeissler fc Masciadril 
2006) This effect makes it difficult to study the reliability of 
the ECMWF data at high altitudes in winter. In this sea- 
son the balloons frequently explode. This happens, highly 
probably, because, due to the low pressure in the high part 
of the atmosphere, the balloons increase their size and they 
finally explode. The elastic material of the balloons is much 
more fragile in winter, when the temperature is much lower 
than in summer in the high part of the atmosphere. 

In local summer the median difference in the wind speed 
never exceeds 1ms -1 . Closest to the ground the difference is 
even smaller, below 1 m s . During local winter the median 
difference in wind speed never exceeds 0.5 m s _1 in the upper 
atmosphere, though the radiosoundings only reach ~10 km 
above the ground. The median difference is larger closest to 
the ground (~ 3 m s _1 ) in winter. 

The median absolute temperature difference in summer 
is below 2 K throughout the whole 20 km. In the proximity 



of the surface this difference is of the order of 1 K, closest 
to the surface even less. During the winter the median dif- 
ference is similar to that of the summer in the high part of 
the atmosphere. In the first 100 m the median difference is 
significantly larger, more than 6 K nearest the surface. 

Figure [5] shows the same outputs of FigfT] but calcu- 
lated for South Pole. We re port just the first 150 m bec ause 
we had already calculated (|Geissler fc M asciadri 200(J) the 
difference ECMWF analyses/radiosoundings above 150 m 
at South Pole. For the wind speed the conclusion is simi- 
lar to what observed at Dome C. The wind speed difference 
remains within 1 m s _1 in the first 150 m in summer. How- 
ever, analyses show a tendency in overestimating (~ 2 m 
s -1 ) the wind speed in the very low levels in winter even if 
this effect is slightly weaker than above Dome C. The me- 
dian absolute temperature difference in summer is similar 
to what observed above Dome C and contained within ~ 
1 K. Near the ground the ECMWF analyses are almost 2 
K warmer than the radiosoundings. This same trend is ob- 
served also in winter. However, in this season, the analyses 
are visibly much warmer (~ 6 K) than the radiosoundings 
near the surface. The statistic uncertainty cr/yN for the 
wind speed is of the order of 0.2 m s^ 1 all along the 20 
km with a maximum of 0.4 m s _1 at 5 km in summer. The 
statistic uncertainty for the absolute temperature is of the 
order of 0.3 K and slightly larger (~ 0.6 K) in the first 30 
m in winter. The precision of the tempearture provided by 
radiosoundings is ~ 0.2 K in the boundary layer and ~ 0.4 
K in the free atmosphere while the precision of the wind 
speecQ is of the order of 0.15 m s _1 . 

Summarizing, we should say that the wind speed is al- 
most well reconstructed by the ECMWF analyses with ex- 
ception of the surface in winter where ECMWF analyses 
show a tendency in overestimating of 2-3 m s~ . For a typ- 
ical wind speed of 3 m s" 1 this corresponds to a large dis- 
crepancy. The absolute temperature, is in general warmer 
from the ECMWF analyses than from the radiosoundings 
near the surface in winter achieving a difference of the or- 
der of ~ 6 K. The wind speed and temperature show sim- 
ilar trends above the two sites with exception of the abso- 
lute temperature in winter. At South Pole, the temperature 
from ECMWF analyses appears warmer on the whole 150 
m while, at Dome C, only in the first 20 m. 

To check if any biases are present in the measurement 
provided by radiosoundings in the first vertical grid point 
(critical region for radiosoundings) we calculated (Table [2| 
the median values of wind speed and absolute temperature 
in the three central winter months obtained from the Auto- 
matic Weather Station (AWS (An automatic weather 
station is an automated version of the traditional 
weather station, to enable measurements from re- 
mote areas) and we compared these outputs with those 
obtained from radiosoundings in the first grid point. Data 
of 2004 are used for South Pole because no AWS data 
from 2006 were available. The median wind speed from ra- 
diosoundings in first vertical grid point is 2.8 ms" 1 and 5.7 



4 More precisley, the velocity uncertainty is defined in the tech- 
nical specification as the standard deviation of the differences in 
twin soundings. 

5 http:/ /amrc. ssec.wisc.edu/realtime. html 
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Figure 1. The median of the difference of the absolute temperature and the wind speed between ECMWF analyses and radiosoundings 
(ECMWF - radiosoundings) for summer (December, January, February) and winter (June, July, August) at Dome C in 2006. 



Summer: |v| Summer: T Winter: |v| Winter: T 




Figure 2. The median of the difference of the absolute temperature and the wind speed between ECMWF analyses and radiosoundings 
(ECMWF - radiosoundings) for summer (December, January, February) and winter (June, July, August) at South Pole in 2006. 
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Figure 3. The gradients of the potential temperature and the wind speed near the surface for Dome C and the South Pole in July 
2006. The dashed lines refer to Dome C and the solid lines refer to the South Pole. Thick lines are ECMWF analyses and thin lines are 
radiosoundings. In the top plots (a and b) the radiosoundings are interpolated with a step of 100 m, in the bottom plots (c and d) the 
step is 10 m. 

Table 2. The median wind speed and absolute temperature for the winter months (JJA) from the lowest level of the radiosoundings 
compared to values obtained with AWS. First and third quartile are also given. 



vo (ms 1 ) 
1 st quartile median 



j rd 



To (K) 

quartile I s * quartile median 3 rd quartile 



Dome C (2006) 1.5 
South Pole (2004) 4.6 



Radiosoundings 

2.8 3.8 
5.7 7.2 



Radiosoundings 
204.4 208.3 213.1 

205.8 209.7 213.6 



Dome C (2006) 2.0 
South Pole (2004) 3.6 



AWS 
3.0 
5.7 



4.5 
7.4 



204.8 
206.7 



AWS 
209.0 
211.4 



213.6 
215.0 



ms -1 respectively at Dome C and South Pole. These values 
match in excellent way with AWS measurements: 3 ms" 1 
and 5.7 ms -1 respectively at Dome C and South Pole in the 
same periods. We conclude therefore that the measurement 
of the first grid point from the radiosoundings are reliable 
because perfectly in agreement with AWS measurements. 
We note that a similar median wind speed of 2.6 ms" 1 has 
been cal culated at Dome C o n a climat ologic scale (1984- 
2003) bv lAristidi et all (|2005l ) in winter. iHudson fc Brandti 
(2005)calculated at the South Pole a median wind speed of 
~ 5 ms" 1 in winter on a climatologic scale (1994-2003). The 
median value calculated in our paper on a time scale of 1 
year have therefore almost no climatologic trends. If we look 
at the absolute tempearture in Table [5] we conclude that the 



radiosoundings are ~ 1 K colder than the AWS at Dome C 
and ~ 2 K at South Pole. 

Going back to Fig[T]and FiglU this conclusion supports 
and confirms the thesis that we are in front of an overesti- 
mate produced by the ECMWF analysis and not to an ar- 
tifact in the measurements (radiosoundings). However, the 
overestimate is slightly smaller than what predicted for the 
temperature. This means an overestimate of 4-5 K instead 
of 6 K. 

FigJT] and FigH therefore indicate that the ECMWF 
analyses accurately describe the state of the free atmosphere 
above Dome C and the South Pole. In the boundary layer, 
the ECMWF data shows a tendency in overestimating the 
wind speed and the absolute temperature, particularly in 
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winter. Results we obtained indicated that ECMWF analy- 
ses should be treated with more caution. 

These conclusions also confirmed the thesis we had ex- 
presse d in the Introduction concerning the lSadibekova et al.l 
(2006) paper who claimed for a good agreement between ra- 
diosoundings and ECMWF analyses even in the surface. 

In spite of the fact that they used ERA-reanalyses (a 
product having a lower resolution than the MARS catalog 
used in our study) the agreement between radiosoundings 
and analyses in their data matches well with our findings 
and predict a good agreement between radiosoundings and 
ECMWF analyses in summer. Our analysis, extended to 
winter, reveals that, in this season, the agreement is far from 
being good and the sharp changes in wind speed and temper- 
ature closest to the surface measured by the radiosoundings 
are not well reconstructed by the ECMWF analyses data. 

To provide the most comprehensive and compact com- 
parison of ECMWF analyses and radiosoundings above 
Dome C and the South Pole in the proximity of the surface 
we prefer to focus now on the two key parameters defining 
the thermodynamic stability, i.e. the gradients of the po- 
tential temperature and of the wind speed. Indeed, only a 
study of the simultaneous systematic effects on both quan- 
tities can tell us if we can use ECMWF analyses to quantify 
the thermodynamic stability in the surface layer. 

Fig. [3] shows the median gradient of the potential tem- 
perature (left) and the median of the square of the gradient 
of the wind speed (right) in the first 150 m with a vertical 
resolution of 100 m (a and b) and of 10 m (c and d). As 
we can expect, the radiosoundings show a sharper gradient 
than the analyses near the surface. We can observe that the 
ECMWF analyses are able to identify that the gradients 
above Dome C are larger than above South Pole. Unfortu- 
nately, a precise quantification is missing and, even in the 
case of the best vertical resolution (cases c and d), the off- 
set produced by analyses with respect to radiosoundings on 
the two parameters (d6/dz and (dv/dz) 2 ) is not compara- 
ble above the two sites (Dome C and the South Pole). This 
implies that the ECMWF analyses do not smooth out the 
potential temperature and wind speed gradients in a similar 
way above the two sites. 

Knowing that the Richardson number depends on the 
ratio of dd/dz and (dv/dz) 2 we conclude that it is pretty 
risky to draw any conclusion on a comparative analyses of 
the Richardson number in the surface layer between differ- 
ent sites calculated with the ECMWF analyses. As a con- 
sequence we can state that we can retrieve a ranking of the 
three sites with respect to the thermal and the dynamic gra- 
dient in an independent way but we cannot retrieve a rank- 
ing of the three sites with respect to the Richardson number 
in the surface layer using ECMWF analyses. We have there- 
fore to limit us to a comparative analysis of the 88 jdz and 
(dv/dz) 2 . In section 5 we will perform the Richardson num- 
ber comparison in the free atmosphere where we showed the 
ECMWF analyses are reliable. 

For a reliable study of the Richardson number in the 
surface layer we need, at present, radiosoundings. A forth- 
coming paper on this topic is in preparation. 



2.1 dO/dz and (dv/dz) 2 at the South Pole, Dome C 
and Dome A 

As a consequence of the conclusions obtained in the pre- 
vious section, we present here the 'thermal' and 'dynamic' 
properties of the surface layer in an independent way. 

The change of the potential temperature with height 
indicates the thermal stability of the atmosphere. A pos- 
itive gradient is defined as stable conditions, the vertical 
displacement of air is suppressed and thus is the production 
of dynamic turbulence. 

The absence of sunlight in the antarctic night and the 
consequent radiative cooling of the snow surface creates 
a strong temperature inversion close to the ground. The 
monthly median of the gradient of the potential tempera- 
ture in the first 150 m for the three sites, calculated with 
the ECMWF analyses, is shown in Fig. [4] From February to 
October all the gradients are positive and indicate a ther- 
mally stable stratified atmosphere near the surface. The in- 
version above Dome A (thick lines) is particularly intense 
when compared to Dome C (thin lines) and to the South 
Pole (dash-dotted lines) during these months. The value of 
the gradient at the lowest level is significantly larger for 
Dome A almost all the year. From March to August there 
is a very sharp change in the slope of the gradient of Dome 
A at around 20 m above the surface. 

During the summer months all three sites have a neutral 
stratification near the surface, i. e. d9/dz m 0. Vertical mo- 
tion of the air is not suppressed but neither is it encouraged. 
A small perturbation can trigger dynamic turbulence. 

The median of the gradient of the square of the wind 
speed in the first 50 m above the ground is shown in Fig. 
[5] The gradient of the wind speed at the lowest level is 
largest at Dome A for every month (except in June when 
it is slightly larger above Dome C). 



3 RADIOSOUNDINGS: THE SURFACE WIND 
SPEED 

Above the summits of the Internal Antarctic Plateau the 
surface winds are expected to be weaker than elsewhere 
on the Plateau. Fig[6] shows the median wind speed near 
the surface measured with radiosoundings at the South 
Pole (dashed lines) and Dome C (solid lines) from April to 
November. We observed that while it is true that the wind 
speed at the very lowest level is weaker at the peak (Dome 
C) than at the slope (South Pole), there is a sharp increase 
in the wind speed above Dome C in the first few tens of 
meters. At the height of 10/20 m, from May to November 
i.e. in the winter, the wind speed is higher above Dome C 
than above South Pole. Above this height the wind speed 
at Dome C is either higher than or very similar to the one 
observed above the South Pole. 

In the center of the winter (June, July and August), 
the wind speed above Dome C reaches the 8 ms -1 at 20 m 
and 9 ms -1 at 30 m. The sharp change of the wind speed in 
the first 10/20 m matches with our expectations of a large 
wind speed gardient. This is indeed a necessary condition 
to justify the pres ence of optical turbulence in the surface 
l|Agabi et al.l 120061 ) in spite of a very stable thermal condi- 
tions (i.e. a positive gradient of the potential tempearture). 
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d9/dz (K/m) d9/dz (K/m) de/dz (K/m) 

Figure 4. The monthly median of the gradient of the potential temperature for 2005, Dome A (thick lines), Dome C (thin lines) and 
the South Pole (dashed lines). 



iTrinquet et al.l (|2008h , in a small sample of radiosound- 
ings in winter, observed a wind speed of 5 ms _1 at 20 m and 
8 ms _1 only at 70 m. Our results, obtained with a complete 
statistical sample in winter, tells us that the Trinquet et 
al. estimate is too optimistic and we should expect a larger 
wind speed at low heights. 



We conclude that in the first 10/20 m the mechanical 
vibrations, that might be derived from the impact of the 
atmospheric flow, flowing at 8-9 ms _1 on a telescope struc- 
ture, are probably more critical above Dome C than above 
the South Pole and should be carefully taken into account 
in the design of astronomical facilities. 



4 ECMWF ANALYSES: THE WEAKEST 
WIND SPEED IN THE TROPOSPHERE 
AND THE LOWER STRATOSPHERE 

Fig. [7] shows the monthly median of the wind speed from 
ground level up to 25 km a.s.l. for Dome A (blue lines), 
Dome C (green lines) and the South Pole (red lines). During 
summer the wind speed is quite weak in the upper part of 
the atmosphere for all three sites. A maximum is observed 
at roughly 8 km a.s.l. From December to April the median 
wind speed is never larger than 15 m s _1 in the whole 20 
km and above all the sites. As the winter approaches the 
wind speed in the upper atmosphere increases. From April 
to November the wind speed increases monotonically above 
10 km a.s.l. and the highest wind speed is more likely to 
be located in the highest part of the atmosphere. However, 
the rate at which the wind speed grows up is not the same 
above the three sites. The smallest increase rate is observed 
at the South Pole and the maximum rate is found above 
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Figure 5. The monthly median of the square of the gradient of the wind speed in the first 50 m for the year 2005. Dome A (thick lines), 
Dome C (thin lines) and the South Pole (dashed lines). 



Dome C. Differences are far from being negligible. In this 
slab the median wind speed is, at Dome C, almost twice 
that of Dome A or the South Pole. 

iGeissler fc Masciadril l|2006l) found in their investigation 
a similar wind speed trend above Dome C from ECMWF 
analyses for the years 2003 and 2004. The wind speed trend 
above Dome C is therefore confirmed in different years. How- 
ever, as announced in the Introduction, the goal in our pa- 
per is the comparison of Dome C with other sites to evaluate 
which are the best for astronomical applications. 

The different wind speed gradient rate observed in Fig. 
[7] above different sites is highly probable related to the syn- 
optic circulation of the Antarctic regions. Indeed, the jet 
stream, that characterizes the vertical wind speed profile of 
mid-latitude sites at the tropopause level, is absent here. On 
the other hand, the Polar Vortex creates stron g high-altitude 
winds surrounding the Polar High in winter ISchwerdtfegerl 
1 19841 ). Wind speed increases monotonically above 10 km. 
The result observed in Fig[7]can be explained with the fact 



that the South Pole is located near the centre of the polar 
vortex and consequently the influence of the polar vortex is 
weak at this site. Dome A and Dome C are situated further 
from the centre of the continent, and thus from the centre 
of the polar vortex. It is to be expected that the wind speed 
is larger above these sites. The farther from the polar high 
(centre of the Polar Vortex) is the site, the greater is the 
wind speed strength above 10 km. 

According to this explanation the wind speed strength 
in the upper atmosphere should be related to the distance 
of the site from the centre of the polar vortex. If we could 
know the exact position of the polar high we would have a 
perfect tool to identify, a priori, the site with the weakest 
wind speed above 10 km in winter. Unfortunately the Polar 
Vortex is not a perfect cone and the centre of the vortex 
is not located at the same coordinates at all heights. To 
verify the sensitivity to this effect we included in our sample 
a further site (Dome F, (77.31 S, 39.7 E), h=38W m) and 
we added the median wind speed in the same Fig(7] In this 
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picture it is evident that South Pole and Dome C have, 
respectively, the weakest and the the greatest wind speed 
above 10 km. The wind speed above Dome A and Dome F 
is mostly comparable and it is more difficult to state which 
site have the weakest wind speed. Dome A is slightly better 
if we consider a statistical analysis of the whole year but 
the difference is not so important. This indicates that the 
polar high fluctuates probably in a region located between 
the South Pole, Dome A and Dome F as indicates the dashed 
region in Fig[8] 



5 THE RICHARDSON NUMBER FOR THE 
SOUTH POLE, DOME C AND DOME A 

The Richardson number is an indicator of the stability of 
the atmosphere: 



Ri 



86/ dz 



(1) 



9 {dv/dz) 2 

where g is the gravitational acceleration (9.8 m s -2 ), 8 
is the potential temperature and v is the wind speed. The 
atmosphere is considered to have a stable stratification when 
the Richardson number is larger than a critical value, typi- 
cally 0.25. If the Richardson number is less than the critical 



value the stratification is classified as unstable. The smaller 
the Richardson number is the higher is the probability of 
triggering of turbulence. 

The comparison of the Richardson number above dif- 
ferent sites allows us to make a relative estimate of which 
site is characterize d by a higher or lower probab ility to trig- 
ger turbulence. In iGeissler fc Masciadril (I2006l )-Fig. 14 the 
median of the inverse of the Richardson number for Dome 
C are shown in different slabs and periods of the year. From 
such results it is possible to retrieve a comparative analysis 
that necessarily is qualitative. In other words we can say if a 
region shows a higher or smaller probability to trigger turbu- 
lence but we do not have a reference to quantify these differ- 
ences nor can we conclude whether these differences are neg- 
ligible or not. In order to provide insights on this question we 
calculated (Fig. O the median of the inverse of the Richard- 
son number (1/Ri) from Dome C (thick solid lines) and from 
Mt. Graham (thin dashed lines), taken as an example of a 
typical mid-latitude site, for each month during the whole 
year 2005. The inverse of the Richardson number is shown 
instead of the Richardson number itself because the inverse 
shows a better dynamic. The smaller 1/Ri is, the more sta- 
ble is the atmosphere. The 1/Ri is smaller above Dome C 
than above Mt Graham almost everywhere. This result is co- 
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Figure 7. The monthly median wind speed for 2005. Dome A (blue line), Dome C (green line), South Pole (red line) and Dome F (black 
line). 



herent with the opti cal turbulence measurements observe d 
above the two site s l|Egner. Masciadri fc McKennal (|200^ . 
lAeabi et all J2006)) and for t his reason it definite l y con - 
firms the method proposed bv lGeissler fc Masciadril (2006) 
to study 1/Ri as a qualitative relative indicator to rank sites 
with respect to the probability to trigger turbulence. Dur- 
ing September and October, when the median wind speed 
is remarkably strong at high altitudes at Dome C (see Fig. 
0, 1/Ri is actually larger than for the mid-latitude site. 
This confirms, from a qualitative and quantitative point of 
view, that the high part of the atmosphere in September and 
October is a region to be monitored carefully because the 
probability to trigger instabilities is comparable and even 
larger than above mid-latitude sites. A strong increase in 
the wind speed is most certainly the cause of the high 1/Ri 
over Dome C at such altitudes. 

The vertical median profiles of 1/Ri for the three 
Antarctic sites (Dome A, Dome C and the South Pole) are 
shown in Fig. 1101 In local summer the profiles from the differ- 



ent sites are similar to each other. The 1/Ri has a maximum 
at ground level and a smaller peak somewhere slightly above 
6 km a.s.l. Above 10 km a.s.l. the atmosphere is very stable 
for all the three sites. From April/May the instability above 
10 km increases. At the end of the winter (September and 
October) the instability in the upper part of the atmosphere 
is even larger than the maximum value observed near 6 km 
for Domes A and C in summer. The instability of Dome C 
is more pronounced than that of Dome A. Above the South 
Pole 1/Ri shows the best conditions (i.e. the most stable) 
over the whole year. 



6 CONCLUSIONS 

In this paper we provide a first comparison of the atmo- 
spheric characteristics above three different sites on the In- 
ternal Antarctic Plateau: Dome C, Dome A and the South 
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Figure 9. The monthly median for 2005 of the inverse of the Richardson number (1/Ri) for Dome C (thick solid lines) and Mt. Graham 
(thin dashed lines). 



Pole. More precisely we try to answer the specific questions 
defined in the introduction. 

(1) The comparison of the ECMWF analyses with the 
radiosoundings shows that the analyses can accurately de- 
scribe the atmosphere above the Internal Antarctic Plateau 
in the whole range from 10 m to 20 km above the surface. 
During no season does the median difference of the wind 
speed exceed 1ms -1 above the first 10 m. The median dif- 
ference of the absolute temperature is within 2 K in the same 
vertical slab. In the surface layer the wind speed discrepancy 
between analyses and radiosoundings is slightly larger (2-3 
m s _1 ) while ECMWF analyses show a tendency to overesti- 
mate the absolute temperature measured by radiosoundings 
in the lowest level in winter (AT ~4-5 K). A statistic analy- 
sis reveals that most of the radiosoundings explode in winter 
at about 10-12 km. This does not allow us to estimate the 
reliability of the ECMWF analyses at these high altitudes. 

(2) We proved that the ECMWF analyses do not pro- 
duce accurate estimates in the surface layers confirming 



what was only assumed by G eissler fc Masciadril |2006). 
This result represents an answer to our original question ex- 
pressed in the Introduction. The ECMWF intrinsically have 
limitations for the characterization of the atmospheric flow 
in this vertical slab. This justify the employment of meso- 
scale models but, at the same time, also tells us that it will 
be fundamental to prove that meso-scale models can do bet- 
ter than General Circualtion Models (GCM) in this vertical 
slab. 

(3) We could conclude that above all the three sites the 
potential temperature in the surface layer is extremely sta- 
ble even if the ECMWF analyses generally underestimate 
its gradient when compared to measurements obtained by 
radiosoundings. Such an effect is particularly evident in win- 
ter. Dome A is by far the site with the steepest gradient of 
potential temperature and wind speed if compared to the 
South Pole and Dome C. 

(4) We proved that the median wind speed in the first 
meters above the ground is weaker at Dome C than at the 
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Figure 10. The monthly median for 2005 of the inverse of the Richardson number (1/Ri) for Dome A (blue lines), Dome C (green lines) 
and the South Pole (red lines). 



South Pole from April to November. However, the wind 
shear in the surface layer at Dome C is much larger than 
at the South Pole achieving at 10-20 m a wind speed of 8-9 
m s _1 in winter. Such a strong wind shear combined with a 
stable stratification of the air in this layer is most likely to 
be the cause of the intense optical turbulence t hat has been 
meas ured in the first tens of meters at Dome C l|Agabi et al.l 
2006). Such a strong wind speed at this height might be a 
source of vibrations produced by the impact of the atmo- 
spheric flow on telescope structures and should therefore be 
taken into account in the conception and design of astro- 
nomic facilities. 

(5) Median monthly values of the inverse of the Richard- 
son number (1/Ri) indicate that the probability to trigger 
instabilities is larger above a mid-latitude site (for which we 
have a reliable characterization of the optical turbulence) 
than above any of the three sites on the Internal Antarctic 
Plateau. Above all the three Antarctic sites 1/Ri is visibly 
smaller than that measured above Mt. Graham (selected as 



representative of mid-latitude sites). This is the first time 
that such a conclusion has be en achieved and definitely pro ve 
that the method presented in lGeissler fc Masciadril (|2006l ) is 
reliable. 

(6) Moreover, our analysis permitted us also a more 
sophisticated discrimination between the quality of the the 
1/Ri above the three sites. Dome A and the South Pole 
show, indeed, more stable conditions than Dome C above 
the first 100 m. This is probably due to the polar vortex 
which, producing an increase of the wind speed in the upper 
atmosphere, also increases the probability to trigger ther- 
modynamic instabilities. 

(7) We showed that it is risky to retrieve estimates of the 
Richardson number in the surface layer (Fig(3} because we 
did not find an equivalent smoothing effect of the ECMWF 
analyses for the gradient of the potential temperature nor 
for the gradient of the wind speed above different sites. 

(8) In the free atmosphere, above the first 10 km, the 
polar vortex induces a monotonic increase of the wind speed 
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Figure 8. Antarctica map. The sites of South Pole, Dome A, 
Dome C and Dome F are labeled with a black point. The dashed 
region indicates the 'position space' of the polar high at different 
heigths as retrieved from the Fig. [7] 



of the three sites in this vertical slab due to their limited 
reliability in this thin atmospheric slab (see Section 5) and 
radiosoundings are available only for Dome C and the South 
Pole. Above Dome A the gradient of the potential temper- 
ature is particularly large in the very near surface layer in- 
dicating conditions of extreme thermal stability that might 
be associated to a strong value of the optical turbulence 
in this vertical range when a thermodynamic instability oc- 
curs (possibly even larger than above Dome C). However 
our study showed that, to predict the thickness of such a 
layer we should need measurements or simulations with at- 
mospheric mesoscale model with a higher spatial resolution 
near the ground that is able to better resolve the evolution 
of the atmospheric flow. This is a part of our forthcoming 
activities. 

In conclusion, at present, the real solid and unique ar- 
gument that makes Dome C preferable to the South Pole 
for astronomical applications is the extreme thinness of the 
optical turbulence surface layer. We expect at Dome A com- 
parable or even larger optical turbulence values with respect 
to Dome C in the surface layer. We cannot conclude if the 
surface layer at Dome A is thinner than that observed above 
Dome C. However our study clearly indicates that Dome C 
is not the best site on the Internal Antarctic Plateau with 
respect to the wind speed in the free atmosphere as well as in 
the surface layer nor is it the site with the most stable con- 
ditions in the free atmosphere. Both Dome A and the South 
Pole show more stable conditions in the free atmosphere. 



in winter that is proportional to the distance of the site to 
the polar high. Dome C therefore shows the largest wind 
speed above 10 km in winter. At Dome C the wind speed 
at 15 km can easily be almost twice that of Dome A and 
even thrice the wind speed at the South Pole in winter. The 
wind speed above the South Pole is the weakest among the 
three sites on the whole 20 km in all seasons. This conclu- 
sion put therefore fundamental warnings for astronomical 
applications. 

(9) This study allowed us to draw a first comprehen- 
sive picture of the atmospheric properties above the Internal 
Antarctic Plateau. In spite of the presence of generally good 
conditions for astronomical applications, Dome C does not 
appear to be the best site with respect to the wind speed, 
in the free atmosphere as well as in the surface layer. Both 
the South Pole and Dome A show a weaker wind speed in 
the free atmosphere. Estimates related to the surface layer 
need to be taken with precaution. ECMWF analyses can- 
not be used to draw definitive conclusions on comparisons 
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